Px + Wy superfluid from s-wave interactions of fermionic cold atoms 
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Two-dimensional {px + ipy) superfluids/superconductors offer a playground for studying intriguing physics 
such as quantum teleportation, non-Abelian statistics, and topological quantum computation. Creating such a 
superfluid in cold fermionic atom optical traps using p-wave Feshbach resonance is turning out to be challenging. 
Here we propose a method to create &px + ipy superfluid directly from an s-wave interaction making use of 
a topological Berry phase, which can be artificially generated. We discuss ways to detect the spontaneous Hall 
mass current, which acts as a diagnostic for the chiral p-wave superfluid. 
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Introduction: In recent years, the physics of the 2D chiral 
p-wave (px + ipy) superfluids has attracted much attention 
because of its nontrivial statistical properties 12[] and potential 
application in topological quantum computation jslHl- The 
chiral superfluid can also act as a testbed for studying the true 
quantum phenomena such as quantum teleportation and vio- 
lation of Bell's inequality 13], which are often masked by 
the many-body effects in a macroscopic system. There has 
been considerable evidence that the symmetry of the super- 
conducting order parameter in strontium ruthenate (Sr2Ru04) 
is spin- triplet +ipy but the observation of exotic prop- 
erties such as quantum half-vortices and non-Abelian statis- 
tics is a serious problem in Sr2Ru04 because of intrinsic spin- 
orbit couphng in the p-wave order parameter With the 
recent observations of p-wave Feshbach resonances in spin- 
polarized ""^K and ^Li atoms JO, .11^] . a px + ipy su- 
perfluid of fermionic cold atoms may also be reahzable in the 
future 112, T3I [l3l- Schemes for observing anyonic statis- 
tics and implementing topological quantum computation us- 
ing vortices in these systems have recently been proposed 
yl |4|]. A potential advantage of the cold atom px + ipy su- 
perfluid is that, being spin-polarized, non-Abelian statistics 
associated with the Majorana mode in the vortex cores should 
be readily observable. 

In spite of its promise, because of the short lifetimes of 
the p-wave paks and molecules in experiments realiz- 
ing a chiral o-wave superfluid from the p-wave Feshbach res- 
onances [SlIslllOllUl seems, at present, challenging. To cir- 
cumvent this problem, in this paper we propose to use the 
much more commonplace, attractive s-wave Feshbach res- 
onances, coupled with an artificially generated topological 
Berry phase 11511 , to create apx + ipy superfluid of fermionic 
cold atoms. The topological Berry phase, in principle, can be 
generated in a variety of ways. Here we consider the Berry 
phase originating from an effective, artificial spin-orbit cou- 
pling of atoms. In ultra-cold atomic gases, the effective spin- 
orbit coupling can be implemented by ha ving the atoms move 
in spatially varying laser fields iflH 17, 18, 3 20]. Since 
the s-wave Feshbach resonances have already been success- 
fully used to create s-wave superfluids, our method offers a 



promising new way to create a topological chiral p-wave su- 
perfluid directly from the s-wave interactions. We also stress 
that, using the methods and ideas described in this paper, it 
should be possible to realize more complex, chiral d-wave 
{dx2^y2 + idxy) superfluid order parameter starting with p- 
wave attractive interactions. 

Once apx + ipy superfluid is realized in experiments, one 
natural and important question is how to observe the chiral- 
ity of the order parameter In the superconducting state of 
Sr2Ru04, it has been done through the observation of a non- 
zero polar Kerr effect which demonstrates macroscopic 
time-reversal symmetry breaking of the px + ipy order pa- 
rameter However, such methods cannot be used to detect 
the neutral superfluid order parameter, because superfluids, 
in contrast to superconductors, do not directly couple to an 
electromagnetic field. Here we show that it is possible to de- 
tect the spontaneous Hall mass current, a clear diagnostic of 
the Px + ipy order parameter, by coupling the neutral atoms 
to effective 'electric fields' generated by optical potentials of 
laser fields. The method to create a px + ipy superfluid di- 
rectly from an s-wave interaction, coupled with the methods 
to observe the spontaneous Hall transverse mass cuiTent, gives 
a complete description of a promising new way to create and 
analyze a chiral p-wave superfluid in fermionic optical traps. 

Px + iPy superfluid from s-wave interactions: We consider 
N fermionic cold atoms confined to a quasi-two dimensional 
(xy plane) trap. The atomic dynamics along the z axis are 
frozen out by optical traps with a high trapping frequency 
112 ill or an optical lattice with high potential depths 1221 . The 
Hamiltonian of the system is 
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where U (r^) is the external magnetic harmonic trap poten- 
tial on the xy plane. We assume that U (r^) is weak and the 
system can be taken as spatially uniform. V^i (r.^ — rj) = 
gS {vi — Vj) is the attractive s-wave interaction (5 < is the 



2 



interaction strength) between the atoms with opposite spins. 
In this BCS regime, molecules are not energetically preferred, 
therefore their number is strongly suppressed and their effects 
are negligible, /i is the chemical potential and is an ef- 
fective Zeeman field for the atoms. — 7 (p^ x ct^) • 
is the Rashba type effective spin-orbit coupling 1IT6I1 . where 
CTi = (erf, af , af) is a vector whose components are the Pauli 
matrices, and 7 is the spin-orbit coupling strength. 

The spin dependent part /locf + in Eq. ([U, can be en- 
gineered in a variety of ways. For instance, ultra-cold atoms 
moving in a 2D spin-dependent hexagonal optical lattice 1,16,1. 
or in a non-Abelian gauge potential created by spatially vary- 
ing laser fields 11191 l20(l . will experience an effective spin-orbit 
coupling. The spin t and | in the effective Hamiltonian ([T]i de- 
note the effective spins. Their definitions depend on the way 
to create the effective spin-orbit coupling and akeady include 
the spatial dependence of the lasers for generating the effec- 
tive spin-orbit coupling ill d Therefore, the corre- 
sponding laser parameters, such as the intensity and the opti- 
cal lattice spacing, do not appear in the effective Hamiltonian 
([U explicitly, although they do affect the parameters 7 and /lo- 
We note here that our proposal to generate the p-wave inter- 
actions does not depend on the specific methods to generate 
the effective spin-orbit coupling. Because of that, we have not 
specified the definitions of the effective spins in this paper. In 
the experiments, h^, 7, and g can be adjusted by varying the 
laser parameters and the Feshbach resonance. 

Since we consider only the short range attractive interac- 
tion, the Fourier transform of the two-body interaction is ap- 
proximately a constant, V (q) = g. Performing second quan- 
tization we obtain. 
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where |u)'s are the single particle eigenfunctions of the one- 
body part of the Hamiltonian given by Eq. ([T]i, and n is the 
band index. (In the presence of optical lattice potentials, and 
in the absence of spin-orbit coupling, \u) is simply the peri- 
odic part of the Bloch wavefunction). Specializing to the BCS 
reduced Hamiltonian and ignoring the residual interactions. 



(3) 



Note that, for slowly varying scattering potentials 
(r^ — Tj) (i.e., the Fourier transformation V {q) is 
nonzero only for small q = |k' — k|), the multiplicative factor 
on the r.h.s. of Eq. (O can be related to the sum of the Berry 
phases. 
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the up and the down spins II23I1 . However, for the scattering 
potentials which are not slowly varying, V (q) takes nonzero 
value even for large |k' — k| (the situation considered in this 
paper), and the extra phases renormalizing the interaction are 
the Pancharatnam geometric phases I,24i1 . In the following. 



we show that with a Rashba type of spin-orbit coupling, one 
can create an effective px + ipy pairing interaction Vk,k' from 
the s-wave interaction. 

In the presence of Rashba spin-orbit coupling, the one-body 
part of the Hamiltonian can be written as. 
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where is the single particle kinetic energy. The two 
bands of the Hamiltonian ^ corresponding to the eigenval- 
ues E± (k) — Eli — II ± yTi^ + T^fc^ are separated by a gap, 
2ho, at k = 0. We choose /i at the middle of the gap such that 
the Fermi surface, the locus in fc-space where E- (k) — /i, lies 
in the lower spin-orbit band only with the Fermi momentum 
denoted by kp- We limit to the regime •ykp ^ ho ^ Ec (5), 
which necessitates the diagonalization of the spin-orbit cou- 
pling energy as the first step, and work with the lower band 
only, where Ec (g) is the energy scale of the BCS pairing gap, 
determined by the s-wave scattering strength. 

Using the appropriate eigenfunction of Eq. it is 

straightforward to show 
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where / (fc, k') = k'k/Nl (fc') Nl (fc) (a + + fc2)(a + 
\/ oP- + fc'2), N^{k) is the normalization constant for the 
lower-band eigenfunction, 9^ is the polar angle in momentum 
space, and a — ho /'j- 

BCS pairing occurs on the Fermi surface. In the physical 
regime ^kp ':s> ho considered in this Letter, we have kp ^ a- 
In this limit, one can show 
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The interaction in the angular momentum channel m is given 
by u„ (k, k') = ^ d/3yk,k'e""'5, where 13 = 9^,- 9^ is 
the angle from k to k'. We therefore have 



ui (k,k')«ff/4. 
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We see that the bare s-wave interaction at to = chan- 
nel is now completely replaced by the p-wave interaction at 
TO = 1 angular momentum channel in the physical regime 
"jkp ^ ho ^ Ec (g). Consequently, the pairing interac- 
tion is renormalized to a separable interaction in the p-wave 
channel leading to a ground state of a 2D chiral p-wave su- 
perfluid with px + ipy symmetry of the order parameter: 
A(p) = Aq {px + ipy) /pf, where Aq is the energy gap in 
the excitation spectrum of the superfluid. 

The physical origin of the renormalization of the interac- 
tion may be understood through the Berry phase effects of 
a Rashba type of spin-orbit coupling. In the presence of a 
Rashba type of spin-orbit coupling, an atom evolving adi- 
abatically in the momentum space accumulates a geometric 
(Berry) phase associated with the adiabatic change of the mo- 
mentum k, in analogy to the Aharanov-Bohm phase acquired 
by an electron moving in the real space in the presence of a 
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magnetic field. Here the corresponding magnetic field in the 
momentum space is the Berry curvature field 
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The effective "magnetic flux" passing through the Fermi disc 
is $B = J -p^f^k ~ TT as kp ^ a, which means a geo- 
metric phase /3<i>B/27r is obtained for the adiabatic moving of 
atoms from k to k' ( /3 = ^k' — 6'k)- This geometric phase 
is the origin of the additional phase factor cxp (— i/3<i>B/7r) 
in the interaction Vk k' (Eq- IS} around the Fermi surface and 
leads to the p-wave pairing at m = 1 channel. Remarkably, if 
originally the bare interaction is in the px + ipy channel, the 
Berry phase renormalizes the interaction to the d-wave chan- 



nel (m = 2). In this way a 2D 



idxy superfluid 



should be realizable, which is very difficult to create using the 
conventional Feshbach resonance approach. 

In experiments, one can choose a suitable attractive interac- 
tion regime (BCS side) so that the pairing gap for the s-wave 
superfluid would be ^ ?i x 200 Hz. The laser parameters for 
generating the effective spin-orbit coupling should be chosen 
so that the Zeeman field /iq is ~ ?i x 1 KHz. For a typical 
Fermi energy Ep ~ ?i x 1 KHz ifTlll . the spin-orbit coupling 
constant should be chosen so that jkp ^ Ti x lOKHz, which 
should be achievable within the current experimental technol- 
ogy [3 3 l^l- With these parameters, we can limit our 
discussion to the lower spin-orbit energy band and create a 
Px + "iPy superfluid from the s-wave attractive interaction us- 
ing the methods described earlier 

Transverse Hall mass current in 2D px + ipy superfluid: 
Neutral atoms in superfluid can interact with laser fields 
through dipole interactions lEsll . The dipole interaction can 
provide an optical potential, whose gradient can be taken as 
an "effective electric field" for the atoms. Here, we study the 
linear response of a chiral px + ipy fermionic superfluid sub- 
ject to such external effective electric fields which act as a 
perturbation. The following two types of external effective 
electric fields will be considered. 

First, we consider an effective electric field, Ey, applied 
along the y direction, and calculate the transverse response 
of the superfluid along the x direction. This transverse re- 
sponse, which gives rise to a spontaneous Hall mass current, 
is a clear diagnostic of the broken time reversal invariance 
and the associated chirality of the px + ipy order parame- 
ter The transverse Hall current changes the sign as the chi- 
rality of the order parameter is reversed. As such, this mass 
current can be used to detect the realization of the chiral p- 
wave superfluid. In experiments, this effective electric field, 
E = — Vy (r), can be realized by applying a perturbation 
potential V {r) — Vq exp (^~y^/2x^) created by a laser beam 
traveling along the x direction, where x is the beam waist of 
the laser. For simplicity, we set the temperature T — and 
neglect finite temperature effects. We also assume that the ex- 
ternal trap potential is very weak, and neglect the effects of 
the spatial inhomogeneity. 



The antisymmetric component of the spontaneous Hall con- 
ductivity, (Txy = —cTyx, for the chiral superfluid can be ob- 
tained from the anomalous chiral response coefficient, which 
leads to (in momentum and frequency domain) [|26tl . 



axy{ci, uj) ~ / 2hd jv 



(9) 



in the low frequencies to <C Aq region, where d is the thick- 
ness of the superfluid along z direction, h is the Plank con- 
stant, vp is the Fermi velocity. Since we consider a time in- 
dependent perturbation potential, Eq. (|9]l can be simplified to 
Oxy ~ -l/2hd. 

Equation (|9]l can be used to calculate the transverse mass 
current jx (r, t) = j J axy{(l, '^)Ey (q, uj) e^'^''^ e^'^* dqduj in- 
duced by the longitudinal effective electric field Ey. For the 
potential Vb cxp (^—y'^/2x'^), this can be simplified as 



Jx (r) = Voyexp (-y^/Sx') /2hdx^ 



(10) 



We see that at the peak of the potential, y — 0, the mass cur- 
rent is zero. The current flows in opposite directions on the 
two sides of the potential, and reaches maximum aty — ±x- 

In the time-of-flight measurements, such a current would 
lead to a velocity of BCS pairs, Vx (r) ~ jx (r) /n (r), where 
n (r) is the density of cooper pairs. The different velocities 
of the atoms on the two sides yield a larger image along the 
X direction compared to the unperturbed case. The enhance- 
ment of the size of the image is determined by the maximum 
transverse velocity, v"^^^ = Vqc^'^-^ /2nxh, which occurs for 
the atoms aty — ±x. Assuming representative values for the 
parameters, x = 20 /im, Vo/h = lOOHz, the 2D cooper pair 
density n — lO^^m^^, we find v^'^^ — 8/im/s. For a time 
of flight 500 ms, the enlargement of the image is about 4 fim, 
which should be observable in experiments. 

In the experiments, one can also detect the velocity distribu- 
tion Vx (y) directly by the two-photon Raman transition llzsll . 
The Raman lasers are focused on a local region and transfer 
atoms in that region to another hyperfine state |3). Then one 
can detect atoms at the state |3) using the time of flight image, 
leading to a determination of their velocity distribution. 

The second type of effective electric field we consider is 
generated by a laser beam propagating along the z direction 
and centered at {x, y) — (0, 0). The optical potential can be 
written as V = Vb exp {^—r'^/2x^^, where r — ^ x^ +2/^- 
The effective electric field is now along the radial direction, 
which leads to a response current 



3t (r) = K)r exp (-rV2x') /2Mx' 



(11) 



along the tangential direction. Such a current forms a close 
loop around the center that corresponds to a rotation of the 
superfluid (see Fig. [T]l. The velocity reaches its maximum at 
r — X and then decreases on both sides. The direction of the 
rotation is determined by the chirality of the superfluid. This 
phenomenon can be observed by measuring the local veloc- 
ity of the atoms using the Raman process. Remarkably, by 
applying a non-rotating laser beam, one can create a rotation 
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FIG. 1 : Generating a rotating 2D + ipy superfluid through a non- 
rotating laser propagating along the z direction. The figure shows 
the tangential current distribution jt (r) for a static laser potential 
V (r) — Vo exp (— r^/2x^). The current has the maximum at r = 
X and decreases on both sides, where x is the laser beam waist. 



In summary, we have proposed a concrete method to gen- 
erate a chiral px + ipy cold atom fermionic superfluid by ex- 
ploiting the well-established s-wave Feshbach resonance and 
the topological Berry phases, thereby circumventing the short 
lifetime issues of p-wave superfluids associated with p-wave 
Feshbach resonance. We have also proposed techniques for 
the direct observation of the chirality of the neutral Px + ipy 
atomic superfluids in optical traps. 

This work is supported by ARO-DARPA. We thank Prof. 
Qian Niu for stimulating discussion. 



of the condensate in a 2D px + ipy superfluid. Note that the 
total angular momentum is conserved in this process. In a 
Px + iPy superfluid, each Cooper pair carries a unit of internal 
angular momentum. In our scheme, this angular momentum 
comes from the effective spin-orbit coupling for the atomic 
motion. The external non-rotating laser potential produces a 
density gradient of the Cooper pairs, leading to the redistri- 
bution of the angular momentum spatially. This redistribution 
of the angular momentum yields the tangential current peak- 
ing at r — X, and is the physical origin of the rotation of the 
condensate in a 2D px + ipy superfluid. 

We emphasize that there is no antisymmetric transverse 
mass current in an s-wave or px-wave superfluid. Thus, the 
above experiments involving the transverse currents can serve 
as clear diagnostic tests for the existence of the chiral super- 
fluid. In the Px + ipy superfluid, the time reversal symmetry 
is broken, which leads to a non-zero Berry phase in the mo- 
mentum space. The non-zero BeiTy phase, absent in the s or 
Pa; -wave superfluids, is the physical origin of the nonzero an- 
tisymmetric transverse mass cuiTent. 

Finally, we note that our proposed method for observing 
the mass current is very general. It does not depend on the 
specific way to generate apx + ipy superfluid. Therefore, if a 
Px + iPy superfluid can be generated using other methods (say, 
using a p-wave Feshbach resonance), our proposed diagnostic 
methods still apply. 
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